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<^ ! ABSTRACT 

Context. Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds 
detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the 
first glimpse into the nature of these sources. The aim of the project is to derive the physical properties of the full cold core population 
I ■ revealed by Planck. 

Aims. We examine the properties of the dust emission within the three fields observed during the SDP. We determine the dust sub- 
millimetre opacity, look for signs of spatial variations in the dust spectral index, and estimate how the apparent variations of the 
parameters could be affected by different sources of uncertainty. 

Methods. We use the Herschel observations where the zero point of the surface brightness scale is set with the help of the Planck 
satellite data. We derive the colour temperature and column density maps of the regions and determine the dust opacity by a com- 
parison with extinction measurements. By simultaneously fitting the colour temperature and the dust spectral index values we look 
for spatial variations in the apparent dust properties. With a simple radiative transfer model we estimate to what extent these can be 
explained by line-of-sight temperature variations, without changes in the dust grain properties. 

Results. The analysis of the dust emission reveals cold and dense clouds that coincide with the Planck sources and confirm those 
detections. The derived dust opacity varies in the range /c(250//m) ~ 0.05-0.2 cm 2 g _1 , higher values being observed preferentially 
in regions of high column density. The average dust spectral index p is ~ 1.9-2.2. There are indications that p increases towards the 
coldest regions. The spectral index decreases strongly near internal heating sources but, according to radiative transfer models, this 
can be explained by the line-of-sight temperature variations without a change in the dust properties. 

Key words. ISM: clouds - Infrared: ISM - Submillimeter: ISM - dust, extinction - Stars: formation - Stars: protostars 

1 . Introduction pendencies, we need to study a large sample of molecular clouds 

■ and the pre-stellar cores within them. The molecular lines are a 
q A major question in star formation studies is how the main fea- central tool in the study of molecular clouds . However, in the 

tures of the process depend on the initial conditions in cold cloud cores many moleC ules have frozen onto dust grains mak- 
molecular clouds. The masses of the formed stars, the star for- ing it difficult t0 carry out quantitative analysis with mo lecular 
mation efficiency, the differences between clustered and isolated lme data alone (e g ? et al rjggg Qne has t0 resort t0 

star formation, and the timescales are all closely linked to the obser vations of the dust component of the clouds, either using 
properties of dense cloud cores. In order to understand these de- the extinction or the dust emission as a tracer of the mass distri- 
bution. 
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* Planck {http://www.esa.int/Planck) is a project of the European The recentl launched Herschel and Planck satellites repre- 
Space Agency - ESA - with instruments provided by two scientific . • -r , . r A • u n- + ± a- u i- 

v j j u t-oa u r 4-' i *l i j sent a significant step forward m sub-millimeter studies by elim- 

consortia funded by ESA member states (in particular the lead coun- . . r * • m r i • n- 

tries: France and Italy) with contributions from NASA (USA), and tele- matm S the atmospheric effects from observations. The satellites 

scope reflectors provided in a collaboration between ESA and a scien- enable sensitive measurements of the dust emission m Galactic 

tific Consortium led and funded by Denmark. interstellar clouds, Planck over the whole sky and Herschel m 

** Herschel is an ESA space observatory with science instruments smaller regions but with higher angular resolution. The satellites 

provided by European-led Principal Investigator consortia and with im- provide the complementarity needed for a study of a full cross 

portant participation from NASA. section of the Galactic cold core population. Because of its high 
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Table 1. The observed fields 



Target RA (J2Q0Q) DEC (J2000) Map size (PACS/SPIRE) 

PCC 249 22 21 17.6 +63 42 25 50'/50' 

PCC288 22 53 31.3 +62 31 44 18730' 

PCC 550 12 25 16.5 -71 46 03 18730' 



Like the spectral index, also the absolute value of the dust ab- 
sorption cross section, k, is expected to be variable. This could 
be caused by a change in the abundance of different dust popula- 
tions, a change in the grain size distribution (see, e.g., Steinacker 
et al.|2010 ), or a change in the composition of individual grains. 
For the sub-millimetre observations, the effects should again be 
strongest in dense and cold regions where the grains acquire ice 
mantles and may more easily coagulate to form larger grains 
with k values up to 3-4 times higher (e.g., Ossenkopf & Henning 
T9941 Krugel & Siebenmorgen [T9941 Stepnik et al. 120031 An 
error in k automatically translates to a similar fractional error in 
the mass estimates. The absolute value of k is not easy to mea- 
sure because it requires an independent column density estimate 
that is equally or more reliable over the Ay range for which sub- 
millimetre observations exist. 

In this paper we study these questions with the help of pho- 
tometric Herschel observations. We start by deriving the colour 
temperature (Sect. 13.11) and column density (Sect 13 21) maps of 
the three SDP fields. This first analysis is done with a constant 
value of the spectral index j3. In Sect. l3.3l we proceed to study the 
evidence for variations in the dust spectral index. With simulta- 
neous fits of colour temperature and J3, we calculate maps of 
the apparent spectral index paying special attention to the error 
sources that could affect the obtained values or the morphology 
of the J3 maps. By comparing the results with the near-infrared 
reddening of background stars, we calculate maps of dust sub- 
millimetre opacity (Sect. 13.41) . To complement this analysis, we 
present in Sect. |4] a simple radiative transfer model for the field 
PCC288 and use the results to estimate the difference between 
the apparent /3 and the real spectral index of the dust grains. After 
a discussion in Sect.[5]we present our final conclusions in Sect. [6] 

2. Observations 

The Planck satellite is performing all sky surveys at nine wave- 
lengths between 350/mi and 1cm (Tauber et al. 2010). With 
the detection methods described in Montier et al. (2010), three 
target fields were selected from the Planck First Light Survey 
(see Table [B for photometric observations with the Herschel 
PACS and SPIRE instruments (for details about Herschel, see 
Pilbratt et al. l20T0l Poglitsch et al. l20T0l and Griffin et al. 120101) . 
The Herschel observations were performed in November and 
December 2009. The field PCC249 was observed in the par- 
allel mode and the fields PCC288 and PCC550 in the norma l 
scan mapping, one instrument at a time (see Juvela et al. 120101) . 
The fields are shown in Fig. [T] The data were reduced with the 
Herschel Interactive Processing Environment (HIPE) v. 2. The 
PACS maps were created using the madmap algorithm. The pro- 
cessing includes some high pass filtering. In order to minimize 
the removal of large scale structure, the length of the filter win- 
dow was kept equal to the scan length. The SPIRE maps are the 
product of direct projection onto the sky and averaging of the 
time ordered data. 

In order to study dust properties over the whole fields we 
need to establish a zero point for the intensity scales. We use the 
Planck maps (data release DX4, available within the Planck core 
teams) at 350 /mi and 550yum and the Improved Reprocessing 
of the IRAS Survey (IRIS, Miville-Deschenes & Lagache l2005l) 
maps at 100/zm. The highest frequency Planck channels have 
been calibrated with FIRAS (Piat et al. 2002) and the IRIS data 
have been calibrated according to DIRBE. When necessary, we 
interpolate these reference data for the Herschel bands using a 
modified black body curve B v v 2 fitted to the nearest two wave- 
length points. 



sensitivity and the coverage of the sub-millimetre wavelengths 
dominated by cold dust emission, the Planck all-sky maps pro- 
vide an unbiased census of cold dust clouds and, in particular, 
of compact pre- stellar cloud cores. The survey of cold compact 
dust clouds is one of the goals in the Galactic science programme 
of Planck. In the Herschel Open Time Key Programme Galactic 
Cold Cores, we use the Herschel PACS and SPIRE instruments 
to map about one hundred target fields that are selected on the 
basis of the Planck survey. This study was begun during the 
Herschel Science Demonstration Phase (SDP) when three fields 
were observed. 

The main features of the compact sources found in the SDP 
fields were reported in Juvela et al. (120101) . In this paper we con- 
centrate on the dust emission at large scales. We are interested in 
the cloud properties, their temperature and column density distri- 
butions. However, we also want to study the dust component it- 
self. These are not unrelated questions because the derived cloud 
masses depend on the assumptions of the dust opacity and spec- 
tral index. Conversely, the variations in dust properties should be 
related to the physical conditions within the clouds. 

The dust spectral index and its temperature dependence have 
received much attention because they provide additional infor- 
mation on the chemical composition, structure, and the size dis- 
tribution of interstellar dust grains (e.g., Ossenkopf & Henning 
IT9941 Krugel & S iebenmorge n [Tggl Mennella et al. 1998; 
Boudet et al. [20051 Meny et al. [20071 Compiegne et al. I20TT1) . 
The observations consistently show an inverse T - /3 relation 
(e.g., Dupac et al. 120031 Hill et al. 120061 Veneziani et al. l20T0l 
and the initial Herschel results including Anderson et al. 2010 
and Rodon et al. 2010). However, the reliability of this relation 
is notoriously difficult to estimate because any noise present in 
the measurements produces a similar anticorrelation (Schwartz 
et al.[T982l Dupac et al. l2003l Shetty et al. l2009aD . Shetty et al. 
d2009bl) studied the importance of observational noise and tem- 
perature variations on the reliability of the derived dust proper- 
ties. In the presence of line-of- sight temperature variations, the 
observed dust spectral index J3 will underestimate the real spec- 
tral index of the dust grains. The estimated colour temperatures 
are not only biased towards the warm regions but, for example 
in the two component models of Shetty et al. ( 2009b), could be 
even higher than any of the real dust temperatures. Similar re- 
sults were obtained by Malinen et al. (2011) who studied syn- 
thetic observations produced by combining MHD modelling and 
radiative transfer calculations. They concluded that the presence 
of protostellar radiation sources can further increase the bias in 
the observed f3 values and can produce an apparent inverse J3-T 
relation that is difficult to separate from any intrinsic relation the 
dust grains may have. Malinen et al. (120111) also noted that when 
the spectral index varies over the wavelength range included in 
the spectral energy density (SED) fit, the derived spectral index 
can rise even above the real fi(A) of the dust anywhere in that 
interval. Because an error in the spectral index is always asso- 
ciated with a compensating error in the temperature, this affects 
the mass estimates and could even bias the derived density pro- 
files of compact objects. 
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Table 2. The uncertainties of the map zero levels 



Target 


Wavelength 
(pim) 


Zero level uncertainty 
MJy sr" 1 


PCC288 


100 


14.3 




160 


10.3 




250 


3.1 




350 


1.2 




500 


0.4 


PCC550 


100 


0.8 1 




160 


1.4 1 




250 


0.43 




350 


0.23 




500 


0.08 


PCC249 


100 


4.4 




160 


2.6 




250 


2.9 




350 


0.7 




500 


0.21 



Zero level calculated using map averages. 



The lOOyum band contains a contribution from very small 
grains (VSG). The effect is approximately similar in the IRIS 
and PACS 100/mi bands and does not affect the calibration 
comparison at that wavelength. However, the presence of VSG 
at lOOyurn could bias the values interpolated to 160//m and 
250yum. Therefore, before the interpolation, we subtract from 
the IRIS lOOyurn data the VSG contribution estimated with the 
Desert et al. (119901) dust model. The radiation field is scaled so 
that the model reproduces the observed temperature of the big 
grains (BG) that dominate the emission wavelengths longer than 
100 yum. The predicted ratio of VSG and BG emission at 100 yum 
is used to correct the IRIS data. Because the ratio depends on 
the estimated BG temperature, the correction has to be done it- 
eratively. Finally, we do a least squares fit to Herschel surface 
brightness (convolved to a resolution of 4.3') vs. Planck/IRIS 
and use the offset to set the zero level of the Herschel maps while 
keeping the original gain calibration. For the PCC550 PACS 
observations the fits remained uncertain because of the small 
amount of surface brightness variation. In this case the zero point 
was estimated using the difference in the average surface bright- 
ness rather than the extrapolation to zero Planck/IRIS surface 
brightness. 

In the comparison between Planck and IRIS, the data should 
be colour corrected because this could affect the quality of the 
least squares fit. The colour corrections were calculated for mod- 
ified black body spectra B v (T)v 2 that were fitted to the observa- 
tions after a preliminary estimation of the zero points. All ob- 
servations, including Planck and IRIS data, were colour cor- 
rected before repeating the least squares fits. The corrections 
can be applied consistently only if the offsets of the Herschel 
and Planck/IRIS maps are already correct and therefore also this 
correction needs to be done iteratively. However, the colour cor- 
rections are less than 10%, and a sufficient accuracy is reached 
already on the second iteration. The subsequent analysis was car- 
ried out using these colour corrected observations. 

The statistical uncertainties of the offsets were estimated 
with the bootstrap method and are listed in Table O The uncer- 
tainties are of the order of 1 MJy sr -1 but smaller at the longest 
wavelengths with the lowest surface brightness. The error es- 
timates are larger for the PACS channels apart from the field 
PCC550 where, as explained above, the quoted offsets are cal- 
culated using average surface brightness values. 



3. Observational results 

In this section we present the colour temperature, column den- 
sity, and spectral index maps. The relation between these and the 
intrinsic dust properties is discussed in the later sections. 

3.1. Colour temperature maps 

The colour temperature maps of the large grain emission were 
estimated using the lOOyurn, 160 yum, 250yum, 350yum, and 
500 yum observations. The colour temperatures were first calcu- 
lated with a constant value of the spectral index, ft = 2.0. The 
validity of this approximation is investigated in Sect. 13.31 The 
use of a constant J3 can be justified even when the assumption is 
not strictly true because it enables more robust estimation of one 
colour temperature that characterizes the fields. As discussed in 
Sect. |H the colour temperature is often an imprecise tracer of the 
physical dust temperatures. 

All data were smoothed to the resolution of the SPIRE 
500 yum map (37") using a Gaussian kernel. As discussed in 
Sect.[2l the determination of the zero points of the surface bright- 
ness scales and the correction for the VSG contribution were 
done iteratively together with the calculation of the colour tem- 
peratures. In the fit the relative weight of the data are set as- 
suming 10% uncertainty for the SPIRE and 20% for the PACS 
measurements. The temperature maps are shown in Fig. [2 

3.2. Column density maps 

The column densities were calculated using the modified black- 
body fits that resulted in the colour temperature maps of Fig. [2J 
Like the spectral index, the dust absorption cross section may 
vary from region to region but the analysis was performed using 
a constant value of /C85o=0.02cm 2 g" 1 (see Juvela et al. 12010b . 
With a v 2 dependence this corresponds to #250=0. 23 cm 2 g" 1 . 
The column density maps are presented in Fig. [3] where we 
also show 9 arcmin diameter circles around the positions of the 
Planck detections. These were used in Juvela et al. (2010| to es- 
timate the masses of the regions. 

The column densities were transformed into mass with the 
distances adopted by Juvela et al. (120101) , 800 pc, 225 pc, and 
800 pc for PCC288, PCC550, and PCC249, respectively. In the 
same order, the total masses of the areas displayed in Fig. [3] are 
-890 M Q , -31 M Q , and -3800 M Q . The masses above the given 
column density contours, and inside the 9 arcmin apertures, are 
listed in Table [3] The statistical errors are small and the mass 
uncertainty is dominated by the uncertainties of the calibration 
(-10%), the dust opacity (up to a factor of two), and the cloud 
distance (typically -30% corresponding to a -50% uncertainty 
in the mass. 

3.3. Dust spectral index 

The 100yum-500yum data were also fitted with modified black- 
body curves B v (T)v^ keeping both the colour temperature T and 
the spectral index (3 as free parameters. We masked out the bor- 
ders of the maps (~ 3') and the northern end of the field PCC288 
wher e the column densi ty map and the map of dust opacity (see 
Sect. 13.41 and Sect. 15.31) show a suspicious feature. The colour 
temperature and spectral index maps are shown in Fig.|4]and rep- 
resentative SEDs in Fig. [51 The median values of J3 are 2.21 , 2.21 , 
and 1.86 for PCC288, PCC550, and PCC249, respectively. The 
same fits without the lOOyurn data are shown in Appendix IB. 2[ 
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Fig. 2. Dust colour temperature maps for the three SDP fields. The calculations assume a fixed value of the spectral index, fi = 2.0. 
The map resolution is ~37" and the contours are drawn one Kelvin apart. 
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Fig. 3. Maps of the hydrogen column density N(R) derived with the colour temperature maps of Fig. [2] The maps have a resolution 
of ~37" and the contours are drawn at levels where logioN(H) equals 21.5, 21.8, and 22.1. The circles show the 9 arcmin apertures 
used in Juvela et al. (120 101) . 



Table 3. The mass estimates of the observed fields or parts of the 
fields that are defined by a column density limit or correspond to 
9 arcmin apertures (regions P, PI, and P2). 



Field 


Region 


Mass (M©) 


PCC288 


full field 


1108 




fogioAf >21.8 


785 




logioN > 22.1 


242 




P 


288 


PCC550 


full field 


36.2 




log l0 N > 21.8 


13.9 




PI 


12.2 




P2 


11.3 


PCC249 


full field 


3802 




log l0 N > 21.8 


1763 




log m N > 22 A 


627 




PI 


252 




P2 


259 



In Fig. Owe show the relations between the derived spectral 
indices and the colour temperatures. For three randomly selected 
pixels, the black contours show where the x 2 value of the SED 
fit has doubled from its minimum value. 

To get the first idea of the possible effect of the calibration 
uncertainty, we compared two cases where we add errors to the 
250 /mi and the 500 fim maps and recalculate the colour tem- 



perature and spectral index maps. The wavelengths were chosen 
because their relative intensities are important in determining the 
spectral index. The map offsets were first changed by the values 
listed in Table [2] and the surface brightness values were scaled 
corresponding to a 10% gain uncertainty. The changes at the two 
wavelengths were made in opposite directions in order to have 
the maximum effect in the derived ft values. By applying the ~ Icr 
shifts both in the gain and in the offset and at both wavelengths 
simultaneously, we should get a conservative estimate of the un- 
certainties. In Fig.O the resulting shifts of the observed (T,j3) 
pairs are shown for the three random pixels. The full maps of the 
maximum and minimum J3 values are shown in the Appendix, 
Sect. ED 



3.4. Dust opacity 

We calculated extinction maps at a 1.5' resolution using stars 
from the Two Micron All Sky Survey (2MASS, Skrutskie et al. 
[20061) and the NICER method (Lombardi & AlvesllED. The 
values of Ay are obtained from NIR colour excess assuming 
an extinction law with Ry=3A. The zero level of the extinction 
was set with the Schlegel et al. (1998) extinction map. The ex- 
tinction estimates enable us to estimate the sub-millimetre dust 
opacity relative to the visual extinction. The extinction is con- 
verted to total ISM mass with the Bohlin et al. (11978b relation, 
N(Ul + H 2 ) ~ 1.87 x 10 21 Ay that is multiplied by 1.4 m H to ac- 
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T (K) T (K) TOO 

Fig. 6. The observed temperature dependence of the spectral index in the three fields. The values are shown for individual pixels 
(green dots). For three selected pixels, the black contours show the limit at which the^ 2 value of the SED fit has doubled from its 
minimum value. For the same pixels the smaller symbols (joined with solid lines) show the extreme values allowed by the calibration 
uncertainty (see text). The error regions for these extreme cases are drawn with dashed and dash-dot lines. 



count for the mass of helium. The conversion is done for easier 
comparison with the literature but the conversion factor itself is 
appropriate mainly for diffuse medium and may vary across the 
fields. 

Figure |7] shows the 250/mi opacity maps derived from 
the SED fits done with j3 = 2.0. The median values in the 
fields PCC288, PCC550, and PCC249 are 0.11, 0.064, and 
0.061 cm 2 g" 1 . The opacities obtained with fixed value of ft = 
2.0 and the free j3 are identical to within 20%. In the cen- 
tral part of PCC288 the values rise close to A:(250yum) -0.22 
while in the northern colder region of PCC249 (corresponding 
to the Planck detection PI) the peak value is even higher, above 
~ 0.3 cm 2 g" 1 . The value used for the mass estimates in Sect. 13.21 
was 0.23 cm 2 g" 1 . This is clearly above the median values but 
roughly consistent with the values associated with the highest 
column densities. However, in the field PCC550 the values re- 
main everywhere below K(250fim) = 0.12 cm 2 g" 1 . Our typical 
opacity of 0.1 cm 2 g" 1 corresponds to r(250yum)=2.3xl0~ 25 /H, 
slightly more than twice the value found in diffuse medium 
(Boulanger et al. [T996l) . 

4. Radiative transfer model of PCC288 

To find out the level at which temperature variations along the 
line-of-sight may bias the observed values of the spectral in- 
dex and the dust temperature, we made a simple model of the 
PCC288 region. The model was not intended as a fit to the 
PCC288 observations but is only used to estimate the general 
observational biases in a setting similar to this field. The model 
was based on the column density map shown in Fig.[3]and repre- 
sents a box with a size of 5.4 pc (for d=800pc). The line-of-sight 
density profile was gaussian with a FWHM that was 0.4 pc to- 
wards the densest part of the cloud and increased up to 1.6 pc 
in regions of low column density. The numbers roughly corre- 
spond to the scales of the central clump and the whole field but 
are basically ad hoc values. In particular, the density distribution 
of the model is quite smooth with a maximum density of only 
3.7xl0 4 cm" 3 . Because also the wavelength dependence of the 
spectral index could bias the observed J3 values (see Malinen et 
al. 120111) , the modelling was done with a modified version of the 
normal Milky Way dust model (Draine 2003 ) where J3 was set 
to a constant value of 2.0 for all wavelengths A > lOO/mi. The 
cloud is illuminated from the outside by an isotropic interstel- 
lar radiation field (ISRF; Mathis [T9831 Three blackbody radi- 
ation sources with a temperature of 2000 K and luminosities of 
10, 20, and 50 Lq were added to represent internal heating. We 
show the results for a model where the column densities and the 



strength of the ISRF were scaled so that the surface brightness 
levels match the observations of PCC288 to ~ 30%. 



Because the proper handling of regions near the internal 
sources requires high spatial resolution, the model was sampled 
onto a hierarchical grid with the smallest cell size correspond- 
ing to the resolution of 1024 3 uniform grid. The dust tempera- 
ture distributions were solved with a radiative transfer program 
that is based on our Monte Carlo code (Juvela & Padoan 120031) 
but has been adapted for use with hierarchical grids (Lunttila 
et al., in preparation). The obtained surface brightness maps 
at 160 yum, 250yum, 350yum, and 500 /mi were analysed to de- 
rive colour temperature and spectral index maps. Noise was not 
added to the synthetic observations and all wavelengths were 
given equal weight in the (J, ft) fitting. The results are shown 
in Fig. [8] where the three internal sources are clearly visible (in 
order of increasing luminosity from West to East). The left hand 
frame shows the difference between the observed colour temper- 
ature and the mass weighted average of the true temperature. In 
the region shown, the dust temperature is overestimated by be- 
tween -0.25 K in the most diffuse areas and ~8 K at the position 
of the strongest source. The median error is ~0.8K. An error 
in the temperature is associated with an opposite change in the 
spectral index. In a fit with a fixed value of ft = 2, the median 
error would be only 0.4 K. At a temperature of T = 15 K, this 
would correspond to a less than 15% underestimation of column 
density. If J3 is kept as a free parameter, the 0.8 K bias leads to 
a ~30% error in column density when the estimate is based on 
k at 160yum. This decreases to only -10% when the estimate is 
based on 500 yum. However, close to the sources the errors re- 
main large, at least a factor of two. 



In Fig. [8j the median value of the observed spectral index is 
/?=1.92 which is only slightly smaller than the actual value of 
the dust model, J3=2.0. The values are more biased in areas of 
high column density and the observed J3 is consistently below 
1.9 in the south. At the position of the sources ft is reduced fur- 
ther to a minimum value of 1.56. The observations provide good 
estimates of the spectral index only in the more diffuse regions 
but also at a certain distance from the sources where the internal 
heating apparently helps to keep the temperature constant along 
the line-of-sight. Although the column density is highest in the 
central clump, the bias is higher in the south, further away from 
the heating sources. 
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Fig. 7. Maps of the dust opacity derived from the comparison of the NIR extinction maps and SED fits with a fixed value of the 
spectral index, J3 = 2.0. The values are given for 250/mi and relative to the total ISM mass. The contours are drawn starting with the 
value of 0.05 cm 2 g" 1 and with a step of 0.02 cm 2 g" 1 (PCC550) or 0.03 cm 2 g" 1 (PCC288 and PCC249). 



5. Discussion 

5.1. Colour temperatures 

The calculation of the colour temperature maps for the full 
Herschel fields confirmed the Planck detections of cold dust 
emission. Even in the field PCC249, with active star formation 
and a high average colour temperature, the Planck detections 
correspond to regions where the colour temperatures, T ~ 15 K, 
are ~2 degrees below the values of the surrounding areas. The 
Planck detections are based on the cold dust signature that re- 
mains when the warm component traced by the lOOyum IRIS 
maps has been subtracted (Montier et al. 2010). This does not 
require a low colour temperature for the total emission. The 
separation of the warm and the possible cold component is 
more uncertain in warm regions where the warm component 
is strong and must be subtracted with higher relative accuracy. 
The PCC249 sources are very bright with peak surface bright- 
ness reaching in both PACS and the 250 /mi SPIRE band sev- 
eral thousand MJysr -1 . The Planck observations were able to 
locate the two nearby cold regions despite the fact that, accord- 
ing to Herschel, they are small and thus beam diluted: the Planck 
beams have a significant contribution also from the warmer sur- 
rounding regions. 

In the analysis performed with a fixed spectral index value of 
J3 = 2.0, the minimum colour temperatures were 13.8 K, 13.1 K, 
15.6 K, and 14.8 K for PCC288, PCC550, and the northern and 
the southern areas in PCC249, respectively. These are a couple 
of degrees higher than the values reported in Juvela et al. (2010 ) 
where temperatures were calculated after the subtraction of the 
local background. In spite of the different nature of the fields, 
the minimum temperatures are rather similar. Of course, even 
the background subtracted observations do not give a true pic- 
ture of the coldest temperatures that inside the cores are likely 
to be below 10 K (e.g., Crapsi et al. [20071 Harju et al. 2008b . 
With radiative transfer modelling of the dust continuum obser- 
vations one could derive an upper limit for the central temper- 
ature but once the temperature is reduced close to 6 K, the cold 
dust will be difficult to detect even with Herschel observations. 
In our simple model of the field PCC288, the minimum phys- 
ical dust temperature was only ~2 K lower than the minimum 
colour temperature. It is very likely that in the real clumps the 
central densities rise higher than in our radiative transfer model 
and therefore also reach lower physical temperatures. 



5.2. Cloud masses and dust opacity 

In Juvela et al. (2010) the masses were estimated only for com- 
pact objects, with a fixed circular aperture and with background 
subtraction. The masses were quoted also for a 9 arcmin diam- 
eter aperture and these can now be compared with the values 
of Table [3] that were obtained without background subtraction. 
In PCC288, the estimate has increased by a factor of two from 
140 M© to 288 M© . In the case of PCC550 the background level 
is small and the change is much smaller. The previous values of 
9.2 M© and 8.1 M© for PI and P2, respectively, have been re- 
placed by 12.2 M© and 1 1.3 M©. In PCC249 the intensity levels 
are again high and also the contrast between the emission in the 
selected apertures and the surrounding reference areas is small. 
As a result, the masses obtained without background subtrac- 
tion are much larger. The previous estimates of 90 M© (area PI) 
and 93 M© (area P2) have increased to 252 M© and 259 M©. 
Differences are, of course, not surprising because the new mea- 
surements are for the whole mass along the line-of-sight while 
the old ones correspond to the column density excess within the 
aperture. 

In PCC249, the new values are in better agreement with mass 
estimates from molecular line observations. Wang et al. (2009) 
obtained from 13 CO(2 -1) observations a value of 1700 M© for 
a region that corresponds to our area PI. The estimate was cal- 
culated for a cloud distance of 1.67 kpc which, with the distance 
of d = 800 pc assumed here, translates to -390 M©. The value 
derived from Herschel data is still lower but already -65% of 
this number. On the other hand, the situation has not changed 
much for PCC550 where the masses obtained from molecular 
line studies (Vilas-Boas et al. 119941 see Juvela et al. 120101) are 
also higher than our values. In PCC550-P2 the molecular line 
data gave a mass of 16.5 M© and the -50% discrepancy is still 
not very serious. However, for the core PI (core Mu5 in Vilas- 
Boas et al. 119941) the CO analysis lead to a value of - 33 M© 
that is three times our estimate. The uncertainty of k can ex- 
plain most of this discrepancy. The values derived in Sect. 13.41 
were all below A:(250yum) = 0.15 cm 2 g" 1 with a median value 
of 0.073 cm 2 g" 1 . With *(250/zm) = 0.085 cm 2 g" 1 the mass es- 
timate of PI would be in perfect agreement with the value ob- 
tained from the line observations. The radiative transfer model 
of PCC288 hinted that the column density estimates may be bi- 
ased. In the case of PCC550 the final answer to this question may 
be provided by the detailed radiative transfer modelling that is 
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currently under way (Ysard et al., in preparation). One can also 
question the accuracy of the CO mass estimates. The fractional 
abundances are not accurately known, the observations of the 
cold regions may be affected by some degree of molecular de- 
pletion (e.g., PCC550), and, conversely, the optical depths of the 
employed 13 CO lines can become high enough so that radiative 
transfer effects can bias the results (e.g., self-absorption affecting 
the observed 13 CO/C 18 line ratios). Considering all the sources 
of uncertainty, the mass estimates obtained from the lines and 
from dust continuum are in good agreement. 

The median values of the dust absorption cross sections 
/c(250yum) obtained from the comparison with NIR extinction are 
0.13, 0.073, and 0.055 cm 2 g" 1 for the fields PCC288, PCC550, 
and PCC249, respectively. The peak values are higher by more 
than a factor of two and, in PCC550 and PCC249, are associ- 
ated with regions of higher column density. This behaviour is 
expected and could be explained by an increase in grain sizes. 
In PCC288, there is some indication that the dust opacity has 
increased in the central clump, close to the column density max- 
imum. The increase towards the North corner is very likely only 
an artifact. The results could change noticeably even within the 
limits set by the accuracy of the zero points of the surface bright- 
ness scale. Nevertheless, the results suggest spatial variations at 
least by a factor of two. Most of our /c(250 mum) values are in 
the range 0.05-0.2 cm 2 g" 1 . Scaled with v 2 , this corresponds to 
k(850 yum)=0. 004-0. 017 cm 2 g" 1 . The Ay values were calculated 
assuming an extinction law with Ry=3.l. This is appropriate for 
most areas but not necessarily for the densest cores. By adopting 
Ry=5.5, the extinction values would decrease by -16%. Taking 
into account a possible decrease in the ratio N(H)/Ay the calcu- 
lated dust opacities could increase further by more than 30% 
(see Evans et al. [20091 Draine [20031 and the online table 
This would further strengthen the case for an increasing sub- 
millimetre dust opacity in the cold and dense cores. 

The opacities can be compared with values found in 
other studies. Using COBE measurements, Boulanger et al. 
(1996) derived an effective dust cross section for high lat- 
itude diffuse medium. At 850 /mi, the obtained values are 
K(850yum)=0.005cm 2 g" 1 in agreement with the low end of our 
measurements. In the dark cloud IC 5146, Kramer et al. ( 2003 ) 
found the values to vary from k^q/kv ~ 1.3 x 10" 5 in warm 
regions (20-30 K) to ~ 5 x 10" 5 in colder regions with T ~ 
12 K. With the conversion used in Sect. 13.41 these correspond 
to -0.003 cm 2 g" 1 and 0.01 1 cm 2 g" 1 , in general agreement with 
our range of values. However, there are also dense and cold 
clouds where similar increase of dust opacity has not been ob- 
served. Nutter et al. ( 2008 ) reported observations of a filament 
in Taurus and concluded that these can be modelled accurately 
without any spatial variations of dust properties. Juvela et al. 
(2009) estimated for a Corona Australis molecular cloud fila- 
ment a k value that was almost consistent with values found in 
diffuse clouds. This in spite of the fact that the region had more 
than 20 m of visual extinction and had a minimum colour temper- 
ature of ~ 1 1 K although, on the other hand, the sub-millimetre 
measurements may have suffered from some spatial filtering. 
The results should be tested, possibly with the help of Herschel 
observations, to confirm whether the dust behaviour is this differ- 
ent in apparently similar clouds or whether the results could be 
explained by the difficulty of reaching an accuracy of measure- 
ments that is needed for a reliable determination of dust opaci- 
ties. 



1 http://www.astro.princeton.edu/ draine/dust/dust.html 



5.3. Spectral indices 

We fitted the observed (T,fi) values of Fig. [4] with two analytical 
formulas,/? = A(T /20K)~ a and/5 = (6 + ojT)~ 1 . The parameters 
of the relations are listed in Table [4] The obtained T(J3) rela- 
tions are steeper than reported in Dupac et al. (2003 ) who fitted 
PRONAOS observations using the latter formula. Our relation is 
closer to the results from the Archeops experiment (Desert et al. 
2008 ) where the observations fit a relation (3 = AT~ a with an ex- 
ponent value of a = 0.66 ± 0.054. The error estimates in Table[4] 
do not reflect the true uncertainty but only the accuracy of the 
curve fitting. One source of uncertainty is the accuracy to which 
the very small grain contribution could be subtracted from the 
lOOyurn data. However, the temperature and spectral index maps 
derived without the lOOyurn band are relatively similar to the re- 
sults shown in Fig. [4] (see Appendix IB. 21) . We will next examine 
the reliability of the detection of the spectral index variations in 
the presence of other error sources. 

Calibration errors affect the (T, p) values mostly in a sys- 
tematic fashion. In Fig. [6j we indicated for three representative 
pixels the shifts that would result from a 10% change in the gain 
calibration combined with a change in the offset that is equal 
to the uncertainties listed in Table [2] The calibration errors tend 
to move all points along a trajectory similar to that caused by 
random noise without much affecting the detection of an inverse 
T - p relation. In the field PCC550, the spectral index reaches 
values p -2.35 at the locations of the two cores. In the high- 
P version (calibration errors causing higher 250 /mi and lower 
500 yum surface brightness) the centre of the PCC550 filament 
would reach p ~ 2.7 while the background would remain at 
P -2.4. Opposite changes in the calibration would decrease the 
spectral index values in the filament to a value ~2, still 0.2-0.3 
units above the background. This shows that the increase of dust 
opacity within the PCC550 filament is robust against calibration 
errors although the uncertainty in the absolute value of p is high, 
close to ~0.4. The effects of a different calibration are similar in 
the other two fields, changing the absolute level of p without a 
clear change in the morphology. 

We used Monte Carlo simulations to quantify the noise bias 
and the evidence for a negative /?- T correlation. We added noise 
to the surface brightness measurement and repeated the fit of the 
P(T) = AT~ a relation.The statistical uncertainty of the surface 
brightness measurements, <xs, was estimated from the difference 
between the observed maps and the SED fits with a free p param- 
eter. These are conservative estimates because the residuals are 
affected also by calibration errors. The calibration uncertainty 
<tq was assumed to consist of a 10% uncertainty in the gain 
calibration and the zero point uncertainties listed in Table [2] 
To see how sensitive the extracted values are on the statistical 
noise alone, we carried out one set of simulations where cr s was 
added to the observed maps. However, we mainly compare the 
original observations with reference simulations with a flat P(T) 
relation. These are based on the SED fits that were made with 
P = 2.0. In the first set of simulations <xs and cr c were added to 
the input maps. In the second set of simulations one quarter of 
each surface brightness map was additionally scaled with a ran- 
dom factor corresponding to a 10% relative error. This simulates 
mapping artifacts, cr M , whose effect is fundamentally different 
from both cr s (affecting each pixel independently) and cr c (af- 
fecting all pixels simultaneously). While cr c can shift all the (r, 
p) points along almost a fixed trajectory, without a noticeable 
change in the best fitting P(T) relation, the mapping artifacts 
will increase the dispersion of the (T,p) values along the noise- 
induced P(T) curve. This could be misinterpreted as a stronger 
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evidence of an inverse relation between p and T. The compari- 
son between the SED fits and the surface brightness maps sug- 
gested the presence of some such artifacts (see Appendix©. In 
PCC249, these were below -5% apart from the northern end of 
the 160 /mi map where, within an area less than 5% of the whole 
map, the residual rises between 10% and 20%. In PCC288 the 
artifacts remain below 10% when the northern corner has been 
masked out. The same applies to PCC550 except for the 160yum 
map where the difference to the SED fit is above 10% over about 
one third of the map and reaches values 20-30% over an area 
less than 10% of the full map. By introducing a 10% error in 
the quarter of the mapped area and in all the frequency bands, 
we are certain to obtain a conservative estimate of the mapping 
artifacts, especially in PCC288 and PCC249. 

Figure [9] shows the distribution of the parameters a in this 
Monte Carlo study, each resulting from the simulation of 500 
sets of synthetic surface brightness maps. For each set of fre- 
quency maps, the colour temperatures and the spectral indices 
were estimated pixel by pixel, and the analytical j3(T) formulas 
were fitted to the results. The distribution S 1 corresponds to the 
observed maps with added noise, <xs. The distributions Rl and 
R2 correspond to simulations with a constant value of ft = 2, 
Rl containing cr s and cr c and R2 additionally the mapping er- 
rors cr M . The vertical lines show the parameter values obtained 
directly from the observations. 

Because of the noise-induced anticorrelation between T and 
P, the a values are positive even when the input data correspond 
to a constant value of/3 = 2.0. For the same reason, the dis- 
tributions S 1 are shifted relative to the observed values and are 
narrow only because the effect is almost identical in all Monte 
Carlo realizations. In the fields PCC288 and PCC249, the ob- 
served a values are higher than most values in the Rl distribu- 
tion. If mapping artifacts can be ignored, this indicates a real 
anticorrelation between temperature and the observed spectral 
index. The hypothesis a > gets probabilities 97% and 99% in 
the case of PCC288 and PCC249, respectively. However, with 
the simulated mapping artifacts the reference distributions R2 is 
much wider and the corresponding probabilities are reduced to 
69% and 66%. The detections are no longer statistically signifi- 
cant but the numbers depend critically on the actual level of cr M . 
We believe that the noise values used in the simulation were con- 
servative and the truth lies somewhere between the cases Rl and 
R2. Therefore, the Monte Carlo study gives some, although weak 
evidence for a non-flat P(T) relation. In the case of PCC550, the 
values of a were higher in the simulations than in observations. 
This is possible if the mapping artifacts conspired to produce a 
low a value or may indicate overestimation of the errors. 

The Monte Carlo study showed that the global analysis of 
(T, p) points does not give a conclusive proof of spectral in- 
dex variations. That analysis is not sensitive to very localized 
P variations while it may be significantly affected by artifacts 
that are more likely to exist close to map edges and far from our 
dense sources. In Fig. ffl many of the spectral index variations 
are clearly correlated with actual cloud structures and, therefore, 
are unlikely to be caused by mapping artifacts. In PCC288 the 



spectral index increases from -2.2 at the boundaries to -2.5 in 
the central clump before decreasing to ~1.5 at the location of the 
strongest warm source at the centre of the map. In PCC249 the 
warm sources are similarly associated with spectral index min- 
ima,/? -1.0 near PCC249-P1 and -1.3 near PCC249-P2, while 
the values in the central map otherwise are between 1.8 and 2.2. 
It is clear that warm sources are associated with significantly 
lower p. The evidence for high p values (»2) towards the cold- 
est clumps is somewhat weaker because it comes from larger 
scales that might be affected by mapping artifacts. However, in 
PCC288 and PCC550, T and p are clearly anticorrelated also in 
the regions with T <15 K. In PCC249 the situation is similar at 
large scales but not in the cores PCC249-P1 and PCC249-P2. 
These are very near to strong radiation sources which may have 
affected the spectral index estimates. 

We conclude that there is strong evidence of spatial varia- 
tions in the apparent dust spectral index and the magnitude of 
these changes is likely to be at least A - 0.2 and can be up 
to A - 1 towards warm, internally heated sources. The simple 
model of the PCC288 field showed that these apparent variations 
do not necessarily indicate an equal change in the dust proper- 
ties. In the model, the spectral index of the grains could be 0.2- 
0.3 units higher than what is observed (Sect. |4]). This applies 
to regions of high column density and without internal heating. 
The effect is opposite to the observed correlation between the 
column density and p and further supports the idea of higher p 
values for the dust grains in the cold cores. This behaviour is 
expected on the basis of laboratory measurements and can be 
explained by the properties of the grain materials at low tem- 
perature (Mennella et al. 1998; Boudet et al. 112005 [ Meny et al. 
2007 ). On the other hand, the low values of p that were observed 
towards the warm sources in PCC288 can be fully explained by 
line-of- sight temperature variations. It will be difficult to deter- 
mine whether also the intrinsic spectral index of the grains is 
lower in such regions (see also Malinen et al. , 1201 It . 

The exact relation between the observed and real p depends 
on the properties of the object in question (density structure, 
strength of the external radiation field, nature of the internal 
sources etc.) and the results of our radiative transfer model can- 
not be generalized very far. We have started detailed modelling 
of the observations of PCC550. The results of that study will give 
us a better idea on how the observed spectral index variations, if 
real, are related to the actual dust properties. This may also shed 
new light on the reliability of the mass estimates derived from 
Herschel observations (Ysard et al., in preparation). 

The possible presence of mapping artifacts was found to be 
the main problem in the evaluation of the observational T~p re- 
lation. In new observations larger map sizes and increased re- 
dundancy would help to reduce such artifacts and make it easier 
to evaluate their impact. The situation may be improved with 
better data reduction methods and, in the future, one could also 
use Planck data to check for low resolution artifacts in the long 
wavelength Herschel data. When combined with far-infrared 
data, the Planck observations themselves are very suitable for the 
study of the spectral index variations and will be used to examine 
the situation also in cold cores and clumps (Planck Collaboration 
l20TTai20TTB . 

6. Conclusions 

We have studied Herschel surface brightness measurements of 
three areas where the Planck satellite data indicated the pres- 
ence of significant amounts of cold dust. By analysing surface 
brightness maps at wavelengths between lOOyurn and 500 yum we 
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map the dust colour temperature, dust opacity, and dust spectral 
index variations over the fields. The results lead to the following 
conclusions: 

- The colour temperature maps confirm the Planck detections 
of cold dust. The analysis of the total intensity leads to min- 
imum colour temperature of 13-15 K in the three fields. 

- The results show an increase in the dust opacity towards 
colder and denser cloud regions. The range of values for 
the dust absorption cross section was estimated to be 
/c(250yum) ~ 0.05-0.2 cm 2 g" 1 , in agreement with previous 
studies. 

- The median values of the observed spectral indices were be- 
tween 1.9 and 2.2. We witnessed an inverse relation between 
the colour temperature and the spectral index. A Monte Carlo 
study indicated but did not conclusively prove a temperature 
dependence of ft beyond what could be produced by obser- 
vational errors. However, because the changes are well cor- 
related with column density, we believe the variations to be 
real. Very low values down to/? -1.0 were observed towards 
warm sources but these are likely to be caused mainly by 
temperature variations rather than changes in the grain char- 
acteristics. 

- With the help of a simple radiative transfer model, we exam- 
ined the relation between the observed and intrinsic spectral 
indices. The results show that the intrinsic ft value of the dust 
grains is underestimated in regions of high column density 
and, therefore, the observations may also underestimate the 
real variation of the spectral index. On the other hand, the flat 
spectrum near internal heating sources can be explained by 
the line-of- sight temperature variations without any change 
in the dust properties. 

Appendix A: Frequency maps 

Fi gures lAJl4 A . 3 1 show the 100/zm, 160/mi, 250 /mi, and 350/mi 
surface brightness maps for the three SDP fields. The maps for 
500 yum are not included but their appearance, apart from the 
lower intensity and lower spatial resolution, is almost identical 
with the 350 /mi maps. 

Appendix B: Residual maps at 160//m-500//m 

The quality of the B v (T)yP fits and the maps themselves can be 
evaluated by examining the residual maps at the five wavelengths 
used in the SED fit. These maps are shown in Figs. IB. 1MB. 3[ 

In PCC288 the average residuals are -4.0, 16.9, 6.7, -6.3, 
and 1.4MJysr _1 in the five bands in the order of increasing 
wavelength. The corresponding values are -0.35, 2.4, 2.4, -1.9, 
and 0.4 MJy sr" 1 in the field PCC550 and -3.8, 17.3, -0.7, -1.6, 
and 0.47 MJy sr -1 in the field PCC249. Allowing some errors in 
the gain, the estimated accuracies of the offsets (see Table [2]) 
appear to have been realistic. 

In PCC288 the largest residuals, ~ 100 MJy sr" 1 , are found 
at the centre of the 160 /mi map. This is associated with the 
brightest point sou rce a nd is still less than 10% of the surface 
brightness. In Fig. IB. 31 the residual images are strongly satu- 
rated and one can see both negative and positive values on the 
two sides of a source. Although the residuals reach several hun- 
dred MJy sr -1 these are only at a ~5% level. The positions of the 
sources themselves agree in all frequency maps so that astromet- 
ric errors should not be a significant factor. Unlike in the case of 
the fast scanning speed, the intermediate scanning speed used in 
our observations should not result in a significant smearing of 



the PACS beam. Nevertheless, the residuals could be associated 
with an imperfect deconvolution of the detector time constants. 

B.1. Maps of the spectral index uncertainty 

The derived spectral index values are affected by observational 
noise, errors in the zero point of the adopted surface brightness 
scale, and errors in the gain calibration. Figure \BA\ shows maps 
that can be used to assess the importance of gain calibration and 
zero point errors. As discussed in Sect. I3.3I we modified the 
250 fim and 500 jim surface brightness maps by adding or sub- 
tracting the zero point uncertainties from Tableland by scaling 
the maps by an amount corresponding to a 10% gain uncertainty. 
The changes at the two wavelengths were made in opposite di- 
rections in order to have the maximum effect in the J3 value. 
By applying Icr shifts in the gain and in the offset and at both 
wavelengths critical for the determination of the spectral index 
we should get conservative limits for the spectral indices. The 
maps of minimum and maximum spectral indices are shown in 
Fig. EH 

B.2. Temperature and spectral index maps without 100 jim 
data 

The subtraction of the VSG emission from the lOOyurn data re- 
liedon the ratio of the VSG and BG grains in the Desert et al. 
(1990) dust model. This also affects the 160yum data but only as 
an uncertainty of the zero point of the intensity scale. To check to 
what extent the colour temperature and the spectral index maps 
are influenced by the lOOyurn data, we recomputed the param- 
eters using the 160//m-500yum data only. The resulting maps 
of t he co lour temperature and the spectral index are shown in 
Fig. 1531 

The maps are morphologically similar to those shown in 
Fig- SI The largest differences are in PCC550 where, without the 
lOOyurn data, the spectral index map closely follows the shape of 
the dense filament. The appearance of SW part of that map again 
suggests the presence of some artifact, the J3 values decreasing 
close to one. Compared to Fig. |H PCC249 exhibits somewhat 
lower values of J3 and a higher range of temperatures. 
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Fig. 1. Three-colour figures of the three SDP fields. The red, 
green, and blue colours correspond to the 250 yum, 160yum, 
and 100 yum intensities (PCC288 and PCC249) or the 350/mi, 
250yum, and 160 yum intensities (PCC550). Only the area cov- 
ered by both the PACS and SPIRE instruments is shown. See 
Appendix lAl for images of the individual bands. 
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Fig. 4. The colour temperature and dust spectral index maps. 
The values have been determined by a fit to the observations 
between 100/zm and 500 yum. The circles indicate the positions 
for which the spectra are shown in Fig.0 
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Fig. 5. Spectral energy distributions for selected positions. The 
measurements correspond to a single pixel at the centre of the 
circles shown in Fig. [4] (after convolution to the resolution of the 
500 yum observations). The temperature and the spectral index 
are quoted in the frames, the first values corresponding to the 
upper SED curve. The dashed lines are the fits with a constant 
spectral index value of J3 = 1. 1 or J3 = 2.3. 
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Fig. 8. Results of the radiative transfer model of the PCC288 
region. The left hand frame shows the difference between the 
colour temperature derived from synthetic observations and the 
real mass averaged dust temperature in the model. The right 
hand frame shows the obtained spectral index (3. The contours 
are drawn at AT = 1 K and AT = 2 K and at £=1.8 and £ =1.9. 
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Fig. 9. Results of the Monte Carlo study where the distribution 
of the values (T 9 J3) is modelled with a function J3(T) = AT~ a . 
The vertical lines show the values of the parameter a estimated 
from the observations. The histograms S 1 shows the distribu- 
tion of the a values for observations with the added noise <xs. 
These are compared to simulations with J3 =2.0 containing the 
noise components <xs and cr c (distribution R 1) or with additional 
simulated mapping artifacts, cr M (R2, see the text for details). 
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Fig. A.l. Surface brightness maps of the field PCC288. The zero 
point of the intensity scale has been set by comparison with the 
IRIS and Planck maps. In the lower frames the dashed line cor- 
responds to the size of the upper frames showing the PACS ob- 
servations. 
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Fig. A.2. Surface brightness maps of the field PCC550. 
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Fig. A.3. Surface brightness maps of the field PCC249. 
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Fig. B.2. Residual maps of the surface brightness in the field 
PCC550 at four wavelengths. 
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Fig.B.l. Residual maps of the surface brightness in the field 

PCC288 at four wavelengths when observations were fitted with Fi §- B3 ' Residual maps of the surface brightness in the field 

the spectral index as a free parameter. The resolution is 1 arc PCC249 at me four wavelengths. 
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Fig. B.4. Estimated uncertainty of the spectral index maps. The 
first column shows the derived values (as in Fig. [4]). The second 
and the third columns show the minimum and maximum values 
of J3 consistent the estimated zero point uncertainties and an error 
of 10% in the gain calibration. The contours are drawn between 
1.4 and 2.9 in steps of 0.3 units (white contours starting with 
value 2.3). 
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Fig.B.5. The colour temperature and dust spectral index 
map obtained by fitting the observations between 160 /mi and 
500 yum. 



